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ABSTRACT 

Aims. We calculate the dust emission expected at 3.43 and 3.53 pm if meteoritic (i.e. hydrogenated) nanodiamonds are responsible 
for most of the far-UV break observed in quasars. 

Methods. We integrate the UV flux that hydrogenated nanodiamonds must absorb to reproduce the far-UV break. Based on laboratory 
spectra of H-terminated diamond surfaces, we analyse the radiative energy budget and derive theoretically the IR emission profiles 
expected for possible C-H surface stretch modes of the diamonds. 

Results. Using as test case a spectrum of 3C298 provided by the Spitzer Observatory, we do not find evidence of these emission 
bands. 

Conclusions. While diamonds without surface adsorbates remain a viable candidate for explaining the far-UV break observed in 
quasars, hydrogenated nanodiamonds appear to be ruled out, as they would give rise to IR emission bands, which have not been 
observed so far. 
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1. Introduction 

The spectral energy distribution (SED) of quasars is composed 
of emission lines superimposed on various continuum emission 
components. The optical to near-infrared domain is reasonably 
well reproduced by a power law. The far-infrared and ultravio- 
let spectral regions are both characterised by a broad continuum 
excess with respect to the optical-IR power law. The UV flux 
excess is referred to as the big blue bump (BBB). According to 
general belief, it corresponds to the thermal signature from a hot 
accretion disk orbiting a supermassive black hole. The extension 
in the extreme UV of the BBB is expected to provide the ionising 
flux that powers most of the emission lines. 

A serious problem with this picture, however, is that the BBB 
observed in quasars appears to be too soft to account for the high 
excitation lines fro m the broad emission line region (BELR) (see 
iKorista et al.ll997h . In effect, a marked continuum drop or steep- 
ening takes place shortward of ^ 1100 A, which we hereafter 
term the far-UV break. This break is observed in the composite 
quasar SED derived by iTelfer et akl (|2002|) , as well as in indi- 
vidual spectra of quasars with 0.7 :$ z ^ 2 dBinette et al.ll2005l 
hereafter BOS). A possibility might be that this break is more 
akin to a localised continuum trough, followed by a marked re- 
covery in the extreme UV, which is the energy region responsible 
for the high-excitation emission lines. Various mec hanisms that 
could generate such a trough are summarised by Binette et al.l 
(l2007h . State-of-the-art calculations of the SED for standard 
geometrically-thin, optically-thick accretion disks do not repro- 
duce the observed break satisfactorily, but these models assume 



a stationary disk with a vertical structure supported only by gas 
and radiation pressure. To our knowledge, no detailed SED cal- 
culations have been carried out assuming a non-stationary disk 
or with an accelerating wind. Absorption by intergalactic H i and 
by intergalactic dust ha s been discarded bvlBinette et al.l (l2003h 
and BOS, respectively. lEastman et al.l ( Il983l) proposed that Hi 
absorption by local clouds that are accelerated up to 0.8c could 
generate a steepening of the transmitted continuum. To repro- 
duce modern data, this model would need to be fine-tuned and 
extended to the extreme UV where the flux recovery is expected. 
More recently, BOS propose that the break could be the result of 
absorption by crystalline dust grains local to the quasars. They 
successfully reproduced the position and detailed shape of the 
UV break in SO quasars, out of a total sample of 61 objects 
from HST-FOS archives whose spectra extended down to at least 
900 A (rest-frame). Their model required two flavours of crys- 
tallites: nanodiamonds without surface adsorbates and nanodia- 
monds with surface adsorbates as found in most measurements 
of meteoritic nanodiamonds from primitive carbonaceous chon- 
drite meteorites; see Mutschke et al. (20 0^ hereaf ter MA04) 
and references therein, and also iJones et al.l (12004 hereafter 
JH04). If the far-UV break is due to dust absorption, we expect 
to observe the infrared (IR) re-emission corresponding to the ab- 
sorbed energy. This is not a straightforward test in the case of 
pure cubic diamonds, sinc e only mu lti-phonon modes are active 
in the IR (JH04; Braatz et al.il2(j0(jt iEdwards|[l985h . However, 
impurities or structural defects can enable one-phonon modes 
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to be come active dAndersen et al.|[r998l : iJones & d'Hendecourll 

l20qol). 

iGuiUois et al] (1 19991) and IVan Kerckhoven et all (l2002i 
hereafter VK02) make a convincing case for the detection of hy- 
drogenated nanodiamonds within circumstellar media for which 
the C-H stretch emission bands at 3.43 /vm and 3.53/vm were 
observecQ. As laboratory measurements of meteoritic nanodia- 
mond sjiavestown^^ 

Il989t iHill et al.lll997l: lAndersen et al.l[T99i iBraatz et al.ll200d 
MA04), we can test their presence by using a method similar to 
VK02. We adopt the spectra of the quasar 3C298 as an example 
test-case, since UV and mid-IR coverage are both available from 
the HST and Spitzer archives. 

2. Calculations of IR emission from hydrogenated 
diamonds 

Our first step consists in estimating the far-UV energy absorbed 
by hydrogenated (meteoritic) nanodiamonds, following the pre- 
scription of BOS. The second step consists in defining a geome- 
try for the dust and computing the IR emission profile using the 
same method as VK02. The derived profiles are finally compared 
to the mid-IR observations of 3C298. 

We assume the quasar continuum source to be isotropic and 
surrounded by a dust screen that lies some distance away from it. 
For the sake of simplicity, the screen can be viewed as a spherical 
shell that covers a fraction Q of the source's sky. The continuum 
source presumably consists of an accretion disk. Although the 
dust emits isotropically, the far-side dust screen is not directly 
observable, due to the opaqueness of the (intervening) accretion 
disk, which accounts for the factor 1/2 in Eq.[T]below. The in- 
trinsic SED over the whole UV domain is modelled as a power 
law of index a, with Fy oc v+" and oc To derive a 

more general formulation of the dust absorption trough, we use 
the concept of equivalent width EWau which we define relative 
to the continuum flux at 1350 A, longward of the break, which is 
a spectral region relatively free of lines (e.g. Fig. 5 of BOS). The 
total energy absorbed by the dust and expected to be re-emitted 
in the IR is given by 

pIR ^ pobs E-w., i-n 
em 2 1350 ^'"^Al , y-) 

where /^°35o '■^^ observed flux at 1350 A (rest-frame) for the 
quasar of interest. The value of EWm is given by the integral 

1300 A 

/ {\-T,)F,dA 
EWf,, = ^52^- , (2) 

1350A 

where T,x - exp(-T,i) is the transmittance of the dust in terms 
of its opacity r^, and F^ the UV flux impinging on the dust 
screen. The opacity can be formulated in terms of the screen 
column A^H times the extinction curve cross-section normalised 
to H. In the model of B05, the dust typically consists of 40% 
of nanodiamonds without surface adsorbates and 60% of hydro- 
genated nanodiamonds. Hence, the opacity due to hydrogenated 
nanodiamonds from meteorites is 0.6A^h o"^' , where the latter 
cross-section corresponds to the extinction curve Al calculated 

' Pure diamonds are extremely inefficient emitters compared to the 
one-phonon mode provided by vibrational modes of carbon bonds 
with H atoms, which populate the surface of meteoritic nanodiamonds. 
Therefore, considering only one-phonon modes for the cooling of the 
grains is considered an acceptable approximation. 



by B05 for hydrogenated nanodiamonds, assuming the small 
grain size regime. To describe the SED, we adopt a spectral in- 
dex of a - -0.45, which is essentially the mean value encoun- 
tered by B05 among their 39 class A quasarfl B05 find that the 
H column that fitted the UV break of their class A quasars was 
1.02 + 0.29 X 10^° cm"^. We obtain a satisfactory fit of the far- 
UV break in 3C298 by assuming a column of 1.2 x 10^" cm"^, 
the value assumed hereafter Integrating Eq.|2l we derive a value 
of EWai of 930 A, which we consider representative of the ma- 
jority of the class A spectra studied by B05. 

To get a handle on the dust covering factor Q, we define, on 
the one hand, the maximalist case where the source is covered in 
all directions by the dust screen, and on the second, a minimalist 
case, in which the dust only covers the source's sky as defined 
by the ionising bicone. The assumed bicone aperture is 45°. The 
two cases can be summarised as follows: 

Case A: O = 1, (3) 
Case B: 0=1- cos(4572) = 0.076. (4) 

With case B, we suppose that the UV radiation that emerges out- 
side the ionising bicone is absorbed by the opaque torus rather 
than by crystallite dust. Since the UV break is ubiquitous in 
quasars, we consider the minimalist case B as a plausibly strict 
minimum for Q.. The difference between the two cases amounts 
to an order of magnitude and represents the level of uncertainty 
regarding geometrical considerations. 

3. The quasar 3C 298 as test-case 

The well-known radio-loud quasar 3C 298 has a redshift of z = 
1 .436 that allows the Spitzer IR Observatory to cover the relevant 
3.5 jum region (rest-frame), where emission from C-H stretch 
modes is expected to occur We retrieved the S14.0 IRS pipeline 
data from the Spitzer archive. The exposure times were 240 s 
with the Short-Low 2 (SL2H-3) slit and 120 s with the Short-Low 
1 (SLl) and Long-Low 2 (LL2) slits. The slit widths are 3'.'6, 
3'.'7, and 10'.'6, respectively. The spectral resolving power is 86 
at 5.2 fim (2.1 fim, rest) in SL2, 71 at 8.5 yum (3.5 yum,rest) in SLl 
and SL3, 86 at 14.6/im (6//m, rest) in LL2, and it increases lin- 
early with wavelength in each order Two cycles were averaged, 
then nods were subtracted to remove the background emission. 
Spectra were extracted in standard apertures (SL2: 7'.'2 at 6 fim, 
SL3: 7'.'2 at 8 fim, SLl: 14'.'4 at 12yum, and LL2: 21'.'7 at 16nm). 
All orders and nods are well-matched in flux, indicating accurate 
pointing on the slits and no significant extended emission. 

The resulting spectrum is shown in Fig.[T] Within the re- 
gion of interest, 3^/im, it agrees within 3% wit h the pub- 
lished (smoothed) version of the same spectrum bv iHaas et al.l 
(I2005IP I However, we note that our IRS fluxes are lower by a fac- 
tor 2.5 compared to those quoted in the NASA/IPAC Database 
(NED), which consisted of ISO and Spitzer MIPS photometric 
measurements. Variability of the underlying non-thermal com- 
ponent, rathe r than of the dust, is the most likely explanation 
(IClearvet al.i r2007). 

The FOS-archived spectrum for this object indicates a flux 
value of F"'^-^Q = 1.0 x 10"'^ erg s"' cm^^A"' (corrected for 
Galactic reddening). Using Eqs.[T]and|2] we derived integrated 
flux values ofF'l = 4.7 x lO^^o and 0.36 x lO^^o W cm^^ for the 



- The near-UV coverage provided by the HST-FOS spectrum of 
3C298 was insufficient to constrain a effectively. 

' The flux scaling in their Fig. 1 should read 10''* rather than the mis- 
printed value of 10'^ (M. Haas, priv. comm.). 
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Fig. 1. Spectram of 3C298 and predicted emission bands based 
on the empirical work of VK02. The thick gray Une shows the 
Spitzer mid-IR spectrum of 3C298 around the 3.5 /im region 
(rest-frame wavelengths). A typical errorbar (with cr ~ 0.02 x 
10-19 W cm ^ jim ') in the interval 3.2-3.8yum is much smaller 
than the gray line's thickness. The inset reveals the Spitzer spec- 
trum over a wider spectral domain. The predicted emission from 
hydrogenated nanodiamonds, assuming cases A and B, is shown 
by dotted lines. Continuous lines (labelled A' and B') show the 
convolution of the VK02 profiles with the spectral resolution of 
the IRS spectrograph (R ^ 60) (see Sect.imi. 

Q. cases A and B, respectively. That much energy should appear 
as two IR bands at 3.43 and 3.53 ^m. 

3. 1 . The 3.43 and 3.53 IR bands from HD 97408 data 

To model the nanodiamond emission, we adopted the 8- 
components fit of VK02 of the two bandfl that were observed 
in the ISO spectrum of the stellar disk of Herbig Ae/Be star 
HD 97408. The method used by VK02 was to decompose the 
observed 3 micron region using Lorentzian profiles and to deter- 
mine the peak frequency, width, and relative strength of the pro- 
files by a least-square fitting routine. The other two stars studied 
by these authors show similar features, but with a lower S/N. 
The profile was normalised to reflect the values expected in 
3C298 and is shown in Fig.[T]for cases A and B. The continuous 
lines show the result of convolving the nanodiamond emission 
profiles with the spectral resolution characterising the data on 
3C298 {A/AA ^ 60). 

As shown in Fig.[T] hydrogenated nanodiamond emission 
should be detected clearly, if present. In the case A scenario, the 
prediction is of very strong emission bands (thin line A' peaks 
at - 14 times above the quasar continuum height) that should 
swamp the IR spectrum at wavelengths near 3.5 //m, which is not 
observed. For case B, where the emission only originates within 
the near-side ionising cone of the quasar, the nanodiamond emis- 

VK02 provide a comprehensive and self-consistent model of the IR 
bands, in which the stellar UV radiation heats up the dust. The profile 
of the strongest component at 3.53 fim allowed these authors to infer 
a nanodiamond temperature of ~ 900 K and a grain size range of 10- 
100 A. 



Fig. 2. Spectrum of 3C298 and predicted emission bands based 
on absorption profiles measured from nanodiamonds of the 
Orgueil meteorite, as measured by JH04. The notation is the 
same as in Fig.[T] 



sion is much less dramatic but still comparable to the height of 
the continuum. The lack of any dust signature on the IR spectrum 
of 3C 298 appears to rule out that hydrogenated nanodiamonds 
similar to the meteoritic type are present in the amount predicted 
by the UV absorption model of B05. 

3.2. The 3.43 and 3.53 fim bands from the laboratory 

Independent laboratory measurements of nanodiamonds from 
distinct meteorites show different peak positions and strengths 
for the features around 3.43 and 3.53 jum (see Table 1 in 
IBraatz et al ] l2000h . It is not clear whether these variations dis- 
close a real physical difference between the diamond populations 
from different carbonaceous chondrites, or if they are the result 
of the laboratory procedures. MA04 used 100 g of the Allende 
meteorite, while JH04 used 1 g of the Orgueil meteorite for their 
otherwise almost identical laboratory measurements. The data 
by MA04 show similar features to those observed by VK02, 
while the data from JH04 do not present the same features. JH04 
compared the Orgueil data with the bands produced by larger 
hydrogenated synthetic diamonds and, in contrast to VK02, 
conclude that the spectra observed from HD 97408 and Elias I 
are better matched by the larger diamonds' profiles. Supposing 
that JH04 is right, instead of using the empirical profiles from 
HD 97408, we should then consider using the profiles measured 
by JH04 (their Fig. 1) corresponding to the nanodiamonds of the 
Orgueil meteorite. In this way, we would ensure that the grain 
size regime remains consistent with the B05 model, which relies 
on a dust model containing meteoritic nanodiamonds of radii 
< 0.6 nm. In Fig.|2]we present the expected IR bands, assuming 
the JH04 profiles inferred from the nanodiamonds of the Orgueil 
meteorite. The areas under the profiles and the notation are the 
same as in Fig.[T] As can be seen by comparing both figures, the 
conclusions reached do not depend on whether we follow the 
interpretation of VK02 or of JH04. 
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4. Discussion 

No emission bands are observed in 3C298 at 3.43 and 3.53;um 
(rest-frame), which we could ascribe to hydrogenated diamond 
C-H stretches, whether we assume the VK02 or the JH04 pro- 
files. Before concluding that meteoritic (i.e. hydrogenated) nan- 
odiamonds are ruled out in quasars, we must contend with the 
following issue. We have assumed that all the energy absorbed 
in the UV wa s to reappear entirely within the 3.43 and 3.53 yum 
bands. Hones & d'Hendecourt (2000) differ on this point, as they 
propose that the nanodiamonds account for many of the PAH 
features as well. However, VK02 argue against this view. On the 
other hand, even if we considered the whole of the 7.7 fim feature 
and the 20 yum plateau observed in HD 97408 to be the results of 
nanodiamond emission, this would still leave 30% of the energy 
to be emitted at 3.43 and 3.53 jum. Even after a reduction by a 
factor of three of the area below the profiles in Figs.[T]and|2l we 
estimate that the inferred absence of the bands would still hold as 
a significant result, at least for intermediate values of f2, between 
cases A and B. 

Our conclusion is that hydrogenated nanodiamonds identi- 
cal to those found in meteorites are not present in the mid- 
IR spectrum of 3C298. An interesting possibility, however, is 
that the surface C-H bonds have been photo-destroyed by the 
quasar's intense UV flux. Such nanodiamonds could not gen- 
erate the bands discussed in this paper. Furthermore, the UV 
properties of the crystallites should mimic those of cubic ter- 
restrial diamonds rather than t hose of meteoritic nanodiamonds 
(with surface C-H bonds). iHaro-Corzo et alj ( l2007h modify 
the original B05 dust model by replacing the meteoritic type 
by an enlarged size distribution of the cubic terrestrial type 
alone. This dust model reproduces the far-UV break of quasars 
and would not generate any C-H stretch emission. To verify 
such a model, we n ow need to consider b ulk impurities of the 
C-N or C-O type dAndersen et al. 19991) or the much less ef- 
ficient multi-phonon mo des (Br aatz et al.ll2000l : lEdwardsll 19851: 
iJones & d'Hendecourll200Ql) . 
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